Anomalous cooling of the parallel velocity in seeded beams 
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We have measured the parallel velocity distribution of a lithium supersonic beam produced by 
seeding lithium in argon. The parallel temperature for lithium is considerably lower than the 
calculated parallel temperature of the argon carrier gas. We have extended the theory of supersonic 
cooling to calculate the parallel temperature of the seeded gas, in the limit of high dilution. The 
theoretical result thus obtained is in good agreement with our observations. 
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I. INTRODUCTION 

We have observed an unexpected cooling effect in an 
argon supersonic beam seeded with a small amount of 
lithium: the lithium parallel temperature is smaller than 
the argon parallel temperature. This effect is surprisingly 
large, the lithium parallel temperature being 3 times 
smaller than the calculated argon parallel temperature, 
whereas the common assumption is that these two par- 
allel temperatures are roughly equal. 

Although surprising, this effect can be explained by 
an extension of the theory giving the terminal parallel 
temperature in supersonic beams of pure gases. This 
theory, first developed by Anderson and Fenn in 1965 
has been completed by Toennies and Winkelmann 
in 1977 and Beijerinck and Verster Q in 1981 (for 
a brief review, see ref. Q)- The expansion of binary 
mixtures was studied by Anderson and coworkers 0, Q , 
as well as by other authors quoted by D.R. Miller in 
his review p. These studies were devoted to the case 
of a heavy species seeded in a light species and the fact 
that the terminal parallel temperatures of the two species 
are not equal was established very early (see figure 8 of 
reference |5(). However, we have not found any previous 
theoretical or experimental evidence of a difference in 
parallel temperatures as large as the one we report here. 

The basic phenomenon is due to the exchange of energy 
between the alkali atoms and the rare gas atoms. Two 
effects are important: i) the alkali-rare gas interaction 
has a longer range than the rare gas-rare gas interactions 
so that the expansion cooling goes on for a longer time 
for the alkali atoms; ii) the mass ratio plays an important 
role in the energy exchange between the two species and, 
as shown below, the parallel temperature of the seeded 
gas (mass 7712) is considerably lower than the parallel 
temperature of the carrier gas (mass mi) when 777,2 <C m\. 

In this paper, we first describe our measurement of the 
parallel temperature of lithium seeded in argon. Then, 
we recall the theory explaining the final temperature in 
supersonic beams and we briefly describe its extension to 
the case of a gas mixture in the limit of a high dilution. 
We thus get differential equations coupling the parallel 



and perpendicular temperatures of the seeded gas to the 
same quantities of the carrier gas. Finally, we compare 
our theoretical result to our experiment as well as to an- 
other evidence of this anomalous cooling effect. 



II. OUR MEASUREMENT 

We have measured the parallel temperature of our 
lithium beam by laser induced fluorescence. The fluo- 
rescence of the 7 Li isotope (92.5% natural abundance) is 
excited by a single frequency laser operating at 671 nm, 
using the 2 Si/ 2 — 2 P3/2 resonance transition. Two laser 
beams are used, one near normal incidence measures the 
transverse velocity distribution, while the other one, with 
an incidence close to 45°, is sensitive to a combination of 
the parallel and perpendicular velocity distributions. We 
use a low laser power density in order to prevent any sat- 
uration broadening. Figure 1 presents a typical spectrum 
of the fluorescence as a function of the laser frequency. 
Each laser beam excites two resonance lines correspond- 
ing to the two hyperfine levels F = 1 and F — 2 of the 
2 Si/2 ground state, with an hyperfine splitting equal to 
803.5 MHz, while the very small hyperfine structure of 
the 2 P3/2 excited state cannot be resolved 0. 

The distance between the lines corresponding to the 
same transition excited by the two laser beams gives 
access to the mean velocity of the beam, while the 
linewidths give access to the widths of the velocity dis- 
tributions. The observed perpendicular velocity distri- 
bution depends on the width of the observed part of the 
lithium beam and our observations correspond to a per- 
pendicular temperature in the 0.4 — 0.7 K range. As ex- 
pected, the parallel velocity distribution is broader than 
the perpendicular velocity distribution (in our excitation 
geometry with a 45° incidence, the Doppler width is only 
0.71 times the width we would get if the laser and atomic 
beams were parallel) . The measured parallel temperature 
is Tj| = 6.6 ±0.6 K. Unfortunately, as we cannot measure 
the parallel temperature of the argon beam, we must de- 
duce its value, 19.0 ± 2.7 K, from the source parameters. 
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III. TERMINAL PARALLEL TEMPERATURE 
IN SUPERSONIC EXPANSIONS 

We consider here only the case of monoatomic gases. 
The starting point is the Boltzmann equation. In the 
frame moving with the hydrodynamic velocity u, the ve- 
locity distribution is assumed to be Maxwellian with dif- 
ferent parallel and perpendicular temperatures Tim and 



F « A(z)(7||i - Tii) with 
A(z) = 16n 1 (z)fi^ 2) (T m )/(15 Uoo ) (5) 

Here, fi^" 1 (T m ) defined by reference || (see equation 
7.4-34) is a thermal average at the mean temperature 

T m = (Tj|i + 2Tj_i)/3 of the collision integral Q\ \ given 
by: 
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mi is the atomic mass. We use local cartesian coordi- 
nates, the x and y axis being in the perpendicular direc- 
tions and the z axis being in the parallel direction. The 
theoretical treatment of Toennies and Winkelmann Q 
(noted below TW) and the one of Beijerinck and Ver- 
ster (noted below BV) are different but equations 
(14), (15b) and (16) of reference are equivalent to equa- 
tions (29a, b) and (B6,7) of reference However, small 
terms of the order of S^ 2 , where the parallel speed ra- 



tio S\\ is defined by S\\ = ^Jmu 1 j(2k B T\\ ) are neglected 
in the BV calculations, which are a good approximation 
only when the terminal value of the parallel speed 
ratio is large. As we are interested in this case here, we 
follow the BV calculations because they are simpler. The 
parallel and perpendicular temperatures are given by: 
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where F is due to the collision-induced energy transfer 
from parallel to perpendicular degrees of freedom: 
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Vi and V2 are the atom velocities before the collision and 
g = |vi — V2I is their relative velocity. The indices 1, 1 re- 
call that the collision considered here involve two atoms 
of species 1. AE is the energy transferred during one 
collision from the parallel degree of freedom to the per- 
pendicular ones; after averaging over the azimuth of the 
final relative velocity around the initial relative velocity, 
AE is given by : 



AE = ^[gl+ 9 l-2gl} [l - cos 2 x ] 



(4) 



where x is the deflection angle. F is well approximated 
by a linear function of (Tim — Tli) : 
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When quantum effects can be neglected, if the interac- 
tion potential is approximated by a 12-6 Lennard- Jones 
potential with a Cq(1, 1) long-range coefficient and if the 
relative kinetic energy is small with respect to the in- 
teraction potential well depth e, the fl\ [ (T) integral is 
given by 0: 

^ 2 i 2) (T) = 2.99(2fc B T/m 1 ) 1/2 (C 6 (l,l)/A: B T) 1 / 3 (7) 

Knowing A(z), which behaves like z _2 T^/ 6 , the differ- 
ential equations can be integrated numerically. The 
parallel temperature tends toward a constant while the 
perpendicular temperature decreases indefinitely, a sit- 
uation very far from thermodynamic equilibrium. The 
perpendicular temperature cools for geometrical reasons 
discussed by reference 0- This calculation relates the 
terminal parallel temperature T^ to the density no and 
temperature To in the source and to the nozzle diameter 
d. The tradition is to give the terminal value of the 
parallel speed ratio: 



= A 
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(8) 



Using SI units, the various values of A and a are A = 
1.412 and a = 0.53 obtained by TW @, A = 1.313 and 
a = 0.545 obtained by BV 0, A = 1.782 and a = 0.495 
obtained also by BV from a fit of experimental val- 
ues for argon expansions. We will use this semi-empirical 
form of equation © to calculate the argon parallel tem- 
perature and we have estimated a 7% error bar on 
from the dispersion of the data set used by BV and by 
comparison with other measurements . 



IV. GENERALIZATION TO THE CASE OF A 
MIXTURE OF TWO MONOATOMIC GASES 

Our calculation is made in the case of a high dilu- 
tion, i.e. when the seeded gas density 77,2 is considerably 
smaller than the carrier gas density n\. We neglect the 
velocity slip i.e. we assume the same hydrodynamic ve- 
locity u for both species. Because <C n±, the expan- 
sion of the carrier gas is not modified by the presence of 
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the seeded gas and the equations written above remain 
valid. For the seeded gas, we write similar equations, 
considering only the collisions with atoms of the carrier 
gas. The calculation of the variation of the parallel and 
perpendicular energies of an atom 2 during a collision 
with an atom 1 is straightforward but too complex to 
be detailed here. In order to linearize the equations, we 
have remarked that the products gQi\(g) oc g 1 ^ 3 vary 
slowly with g and we have treated them as independent 
of g. The values of these products are then chosen to 
be coherent with the treatment of the pure gas case de- 
scribed in part III. We thus get the differential equations 
relating Tp and T±2'. 
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with T\\ av = /3Tj|! + aTp, T± av = /3T ±l + aT± 2 , a = 
mi/ (mi + m^) and (3 — mijirnx + m.2). p s is the ratio 
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of collision integrals: 
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collisions, i.e. if the ratio p s is large enough. We have 
plotted in figurc[21the terminal value of the ratio Tp/Tp 
for various values of the ratio p s and for different values 
of the mass ratio m.2 /to i. 



V. COMPARISON WITH EXPERIMENTAL 
RESULTS 

We may now compare our theoretical and experimen- 
tal results for lithium seeded in argon. The argon (from 
Air Liquide stated purity 99.999%) is further purified by 
a cartridge (also from Air Liquide). The lithium pressure 
is fixed by the temperature of the back part of the oven 
at 0.55 mbar. In our experiment with po = 333 millibar, 
To = 943 K and a nozzle diameter d = 200 /im, the par- 
allel temperature of lithium is T|| 2o o = 6.6 ± 0.6 K. Using 
equation JSJ, we deduce the argon parallel temperature 
from our source conditions, Tji loo = 19.0±2.7 K. We thus 
get a parallel temperature ratio Tpoo/Tjiioo = 0.35±0.08. 
This value is in good agreement with our theoretical re- 
sult T||2oo/Tj| loo = 0.31 obtained with p s — 2.55 deduced 
from the Cr values of argon-argon and lithium-argon in- 
teraction [lOj . 

A similar beam with sodium seeded in argon was built 
by D. Pritchard and coworkers 0,0] with a source tem- 
perature near 1000 K, an argon pressure up to po = 3 bar 
and a nozzle diameter d = 70 pm. The measured velocity 
u = 1040 ± 2 m/s corresponds to T = 1039 K and the 
rms velocity width, deduced from the atomic diffraction 
pattern, was found equal to 3.7 ± 0.4% , corresponding 



Po is the ratio of the I = 1 and I = 2 angle-averaged to Tu 2oo = 4.1 ± 0.9 K. From the calculated argon par- 



cross-sections p — Q12/Q1 2- Using the low energy de- 
flection function for a 12-6 Lennard- Jones potential (see 
appendix A of (3), we have calculated p a — 1.32. 

Following BV, we have introduced reduced tempera- 
tures and a reduced z-coordinate. Then, the numerical 
integration of the coupled equations gives the tempera- 
tures represented in figure [3 At the end of the expan- 
sion, the parallel and perpendicular temperatures of the 
seeded gas are lower than the similar quantities for the 
carrier gas. It is difficult to have a very simple explana- 
tion of this effect, because of the complexity of equations 
(|9I10I) . It is nevertheless clear that the parallel temper- 
ature of species 2 remains coupled to the perpendicular 
temperatures when the parallel temperature of species 1 
is no more coupled. 

The parallel temperature ratio Tp/TiM can be substan- 
tially lower than 1. The minimum value for temperature 
Tp is reached when the right-hand side of equation © 
vanishes. Neglecting the perpendicular temperatures, we 
thus get: 
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This limiting value is reached if the number of 1 — 2 
collisions is sufficiently larger than the number of 1 — 1 



allel temperature Tp = 7.7 ± 1.1 K corresponding to 
the largest pressure, we get the parallel temperature ra- 
tio Tpoo/Tpoo = 0.53 ± 0.19, in good agreement with 
our theoretical result is Tpoo/Tpoo = 0.61, using the Cg 
values from reference 10] (p s = 1.67). 



VI. CONCLUSION 

In this letter, we have described a measurement of 
the parallel velocity distribution of a supersonic beam 
of lithium seeded in argon, with a high dilution. The 
measured lithium parallel temperature 6.6 ±0.6 K is con- 
siderably lower than the calculated parallel temperature 
of the argon beam near 19.0 ± 2.7 K. 

We have extended the theory of the velocity distribu- 
tion in supersonic beams to a mixture of monoatomic 
gases, in the high dilution limit. We have established 
differential equations describing the parallel and perpen- 
dicular temperatures of the seeded gas; these equations 
explain the observed effect as a consequence of several 
circumstances: role of the mass ratio in energy exchange 
during collisions of the two species, longer range of the 
lithium-argon interaction as compared to argon-argon in- 
teraction. We have deduced the lowest possible value of 
the temperature ratio Tpoo/Ipoo- This ratio may be 
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quite low when a light species is seeded in a heavy species 
with m 2 <C mi. 

We have compared our theory with experimental re- 
sults corresponding to lithium seeded in argon (our ex- 
periment) and sodium seeded in argon (an experiment 
done by D. Pritchard and coworkers). In both cases the 
experimental and theoretical values of the ratio of the 
parallel temperatures are in good agreement, but the er- 
ror bars on the experimental values are large. A simul- 
taneous measurement of the parallel temperatures of the 
carrier and seeded gases would be useful to reduce these 
error bars, thus providing a better test of the approxi- 



mations done in our theoretical analysis. 
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FIG. 1: Laser induced fluorescence signal as a function of the 
laser frequency: the dots represent the experimental data and 
the full curves are the Gaussian fits to each line. The Fabry- 
Perot used for calibration has a free spectral range equal to 
251.4 ± 0.5 MHz. The fluorescence peaks are labeled by the 
ground state F value (F = 1 or 2) and by a letter correspond- 
ing to the laser beam: A for the beam near 45° incidence, B 
for the beam near normal incidence. 



FIG. 2: The reduced temperatures r = H" 9/11 T/T are 
plotted as a function of the reduced z-coordinate £ = 
2.48S 12/11 z/d with S = 0.813n d (Ce/k B T ) 1/3 (to give an 
example, for our argon beam seeded with lithium, S = 32.4). 
This plot, which is identical to figure 5 of Beijerinck and Ver- 
ster [3] for Tin and Txi, shows the anomalous cooling effect 
for species 2. We have used p = 1.32 and p a = 2.55 corre- 
sponding to lithium seeded in argon. 
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FIG. 3: The calculated parallel temperature ratio Tj| 2 oo/T|| loo 
is plotted as a function of the mass ratio milmi (spanning the 
range from hydrogen to rubidium seeded in argon) for several 
values of the ratio p a , with p fixed, p = 1.32. From bottom 
to top, p s = oo, ps = 2.55 (lithium case), p s = 1.67 (sodium 
case). The experimental points corresponding to lithium and 
sodium seeded in argon are represented by stars. 
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